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Abstract 

A brief introduction is given to the subject of brane induced gravity. The 
5D example is discussed in detail. The 4D laws of gravity are obtained on 
a brane embedded in an infinite volume extra space, where the problem of 
stabilization of the volume modulus is absent. The theory has two classi- 
cally disjoint branches of solutions ~ the conventional and self-accelerated one. 
The conventional branch gives a perturbatively stable model of a metastable 
graviton, with potentially testable predictions within the Solar system. The 
self-accelerated branch, on the other hand, provides an existence proof for an 
idea that the accelerated expansion of the Universe could be due to modified 
gravity. The issue of perturbative stability of the self-accelerated branch is 
obscured by a breakdown of the conventional perturbative expansion. How- 
ever, a certain exact non-perturbative solution found in hep-t h /0612016 ex- 
hibits a net negative gravitational mass, while this mass is positive on the 
conventional branch. This suggest that the self-accelerated solution must be 
non-perturbatively unstable. A proposal to overcome this problem in an ex- 
tension of the original model, that also allows for the quantum gravity scale 
to be unrestricted, is briefly discussed. 



1 Introduction 



The discovery of Refs. |1] that that the present-day expansion of the Universe 
is accelerating, has been confirmed by a number of subsequent efforts. One way 
to parametrize the accelerated expansion is to postulate the existence of a "dark 
energy" component in the Einstein equation: 

= SttG^vIT^T"" + ^r'^ , (1) 

were G^i, = R^^ — is the 4D Einstein tensor of the metric g^y{x), and 

TJ^^^^^^ and Tj^j^"^^ energy (jg^ote the stress-tensors for matter (including dark matter) 
and dark energy, respectively; the latter has to have a negative enough pressure to 
account for the observations. 

On the other hand, one can consider another logical possibility that the accel- 
erated expansion is due to modified General Relativity (GR). Schematically, the 
modified Einstein equation could be written as: 

G^. - IC,,ig, m,) = StiGnTP''^' , (2) 

where }C^u{g,fnc) denotes a tensor that could depend on a metric g, its derivatives, 
as well as on other fields not present in GR. Moreover, )C contains a dimensionful 
constant rric ~ ifo ~ 10^^^ GeV, which sets a distance/time scale = m^^ at 
which the Newtonian potential obtained from ([2]) significantly deviates from the 
conventional one. One concrete example of lC^v{gi^c) is given by the DGP model 
[2j, and will be discussed in detail in the next section. 

At a first sight, there does not seem to be a difference between ([2]) and ([1]), as 
the new term on the l.h.s. of ([2]) could be transfered to its r.h.s. and regarded as 
the "dark energy" component, similar to the one present in ([T]). 

In reality, however, the difference between ([2]) and ([1]) is significant. It is typically 
implied in ([1]) that the dark energy is either due to cosmological constant or a light 
scalar field. This field forms an independent sector of the theory, the dynamics 
of which is not restricted by severe constraints that general covariance imposes on 
tensor fields. In contrast with this, the new term in ([2]) contains metric itself in 
a nontrivial way. As such it is highly restricted by general covariance. This gives 
rise to both the theoretical and observational differences between models of "dark 
energy" and modified gravity. 

It is worth emphasizing that a primary theoretical motivation for the models of 
modified gravity is to evade the S.Weinberg's no-go theorem on the "old" cosmo- 
logical constant problem (CCP) [3], see Ref. [1| for a summary of these discussions. 
This attractive possibility still exists, in principle, in models of brane induced gravity 
with the number of space-time dimensions D > & O [6l [7] , however, many aspects 
of those models are not well-understood, and we won't be discussing them here. In- 
stead we concentrate on the 5D brane induced gravity [2] . The latter does not offer 
a solution to the old CCP, however, it can be used as an example for understanding 
of a new dynamics introduced by modified gravity. 
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2 Brane induced gravity 



An explicit example of the modified Einstein equation ([2]) is provided by the DGP 
model |2] . All the known 4D interactions, except gravity, are thought to be confined 
to a brane that is embedded in a 5D infinite-volume (uncompactified) empty space 
where only gravity propagates. In this setup, the additional term on the l.h.s. of 
([2]) is provided by the 4D extrinsic curvature terms of the brane 

Gf,y - mc{K^^ - g^uK) = 8nGNT^^{x) . (3) 

Here, (7^,^ = R^^, — is the 4D Einstein tensor for the metric that depends on 

both 4D coordinates x^, and the fifth coordinate y, g^u{,x,y); K = g^^^K^^u, is the 
trace of the extrinsic curvature tensor 

K^,, = ^ {dyg^u - V,N, - V.N,) , (4) 

and is a 4D covariant derivative w.r.t. the metric g,y{x,y). We introduced the 
lapse scalar field N , and the shift vector field A*"^ [5] : 

g,, = N, = g,,N- , g,, = + g,, ^ . (5) 

In the DGP model, equation ^ is accompanied by the {/xz^} equation in the bulk, 
and the {/u5} and {55} equations which respectively read 

Gg)=0 for y^O, (6) 
V^K,, = VuK, (7) 
R = K''- K^.K^'" . (8) 

Here G^fj denotes the 5D Einstein tensor for the 5D metric gAB{x, y) (A,B=0,1,2,3,5), 
and g,^{x,y) is its 4D part. Note that the {/i5} and {55} equations, ([7]) and ([HD, 
should be satisfied in the bulk, ?/ 7^ , as well as on the brane, y = 0. 

The above set of equations can be derived by the variational principle using the 
action [2] 



S = I d!^x ^/glR{g4) H — ^ / d'^x dy ^/glTZ{g^) + surface terms. 



(9) 



where g^^y = g,u{x, 0), and g^ refers to the full 5D metric; R and 71 are the four- 
dimensional and five- dimensional Ricci scalars respectively, and M^, stands for the 
fundamental gravitational scale of the bulk theory. The brane is located at y = 
and Z2 symmetry across the brane is imposed. The boundary Gibbons-Hawking 
term should be taken into account to warrant the correct Einstein equations in the 
bulk. The matter fields, that are also omitted here for simplicity, are assumed to be 
localized on the brane. 
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Another form of the action (Q can be given in the ADM formahsm 



S = ^ j d'x^,Rig,)+^ j d'xdy^N{R + K^-K%). (10) 



The above action gives rise to the equations of motion fl3ll6] |8|). where 

m = ^ (11) 

In order for nic to be of order Hq, we need that the bulk quantum gravity scale to 
be M^, ~ 100 MeV . Such a low value of the quantum gravity scale is in no conflict 
with observations, however, it presents an obstacle for a string theory realization 
of the above model. In the next section, we will show how a mild extension of the 
model can relax the constraint on M*, in which case it could take an arbitrary value 
below the 4D Planck mass |9j. 



3 Role of non-linearities 

To better understand the linearized theory of ([9]), it is instructive to consider an 
analogous scalar model [2] 



t/^x(5^0(x, 0))^ - ^ y d^xdy{dA(t>{x, y)Y , (12) 

where the dimensionless scalar field is to mimic the 5D graviton of the previous sec- 
tion. We impose the Z2 symmetry across the y = boundary, and add to the above 
action the coupling of to a source J, also locahzed on the brane, / d'^xdyS{y)J(f). 
To obtain the junction condition, we integrate the equation of motion obtained from 
(fT2|) w.r.t. y, from 0~ to 0"*". The resulting equation written ai y = reads: 

- dl^\y=o - m,dy(P\y=o+ = J / M^, . (is) 

The bulk equation is nothing but a 5D free scalar equation 

d^dA(l) = for y^O. (14) 

It is straightforward to find a solution to the bulk and junction equations. For this 
we turn to the momentum space w.r.t. the 4D coordinates while keeping the fifth 
dimension intact. The (decaying in the bulk) solution reads: 



where the sign "tilde" denotes the Fourier transformed quantities, and we introduced 
the Euclidean 4- momentum square as follows: = p^p^ = —p^ + Pi +P2 +P3 = 
pI + Pi + P2 + pIj and set the notation p = 
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The solution (fT5|) exhibits a number of interesting properties: 

(1) At short distances (i.e., the large momenta p 3> rric) it reduces to a 4D 
solution with the inverse square law, for the force mediated by this scalar. At 
large distances (i.e., small momenta p -C rric) it turns into a 5D solution with the 

force law. 

(2) The expression f|T5|) appears to have two poles, one at = and another one 
at p = = —rric. However, the residue of the former pole is zero, consistent with 
the expectation that there is no normalizable 4D massless state in the spectrum. 

The second pole appears on a non-physical Riemann sheet of the complex p^- 
plane (see details in Ref. [1]), and describes a resonance-like state. The branch-cut 
on the complex p^-plane is due to the term in (fT5l) . and reflects the presence 
of a continuum of the Kaluza-Klein (KK) states. 

(3) Two localized sources on a brane exchange one 5D scalar state the propagator 
of which can be read off f|T5|) . From the 4D perspective what is being exchanged 
is an infinite number of KK states, the 4D couplings of which are suppressed as 
compared to the ordinary 5D theory by the following factor [ 10] : 

(16) 



1 + (m? /ml) ' 

where m denotes the mass of a given KK state. Therefore, the heavier the KK state 
the more its suppression. That is why at distances ^ = one recovers 4D 
Newtonian potential. 

A similar calculation can be performed for the linearized gravitational theory [2] . 
The object of the primary interests here is the one-graviton exchange amplitude be- 
tween two sources. The expression for the resulting amplitude is gauge independent 
and reads as follows: 

^l-graviton(P, V) = o exp(-p|?/|) , (17) 

where 

As in the scalar example, the pole at = has zero residue, and the second pole 
in f|T7j) is on a non-physical Riemann sheet. Therefore, the amplitude vAi-graviton 
describes propagation of a metastable state with the lifetime ~ "m."^, which decays 
into a continuum of the KK modes. 

However, there is a crucial difference from the scalar case. This has to do with 
the numerator of (fT7|l . In the limit nic the numerator does not reduce to the 
analogous expression in GR. The latter takes the form: 

SitgJtI - It-t] . (19) 



4 



The difference between (fT8|) and (fT9|) is due to the fact that a 5D graviton (or a mas- 
sive graviton for that matter) propagates 5 on-shell degrees of freedom (hehcity-2, 
helicity-1, and hehcity-0), while the GR graviton propagates only 2 on-shell degrees 
of freedom (helicity-2 state). The helicity-1 state of the 5D graviton does not con- 
tribute to the one-graviton exchange amplitude f|T8|) because of the conservation of 
the stress-tensor. However, the helicity-0 state does contribute and gives rise to the 
finite difference between f lTS]) and ffTI?]) . Observationally ffTSl) is excluded! This is 
the essence of the van Dam-Veltman-Zakharov discontinuity (vDVZ) [TT]. 

Fortunately, the vDVZ argument fails for observationally interesting sources. 
This is due to nonlinear interactions [121 [131, easily be understood as fol- 

lows fT3]: the longitudinal part of the graviton propagator in DGP contains terms 
proportional to 



rricP 

This term does not contribute to the amplitude ^i_graviton because of conservation of 
the stress-tensor, however, it does contribute already in the first nonlinear correction 
(since the stress-tensor is only covariantly conserved in the non-linear theory). Due 
to the singularity in rric in (!20|) . the perturbation theory breaks down precociously. 
However, this breakdown is an artifact of an ill-defined perturbative expansion - the 
known exact solutions of the model have no trace of breaking [13] . The perturbative 
expansion in powers of gets "contaminated" by another dimensionful parameter 
l/rric, and this leads to its breakdown. 

Under the circumstances, one could either adopt a different type of expansion, 
e.g., an expansion in the small parameter rric [T3l ITi] . or look at exact solutions 
0. Both of these programs have been carried out to a certain extent. Here we 
summarize the main results. 

The model has one adjustable parameter - the distance/time scale Vc = fn~^. 
Distributions of matter and radiation which are homogeneous and isotropic at scales 
> Tc exhibit the following properties: for distance/time scales <C the solutions 
approximate General Relativity (GR) to a high degree of accuracy, while for scales > 
Tc they dramatically differ [T^l EDI [13] • Since ~ i^o ~ 10~^^GeV, the deviations 
from GR could lead to observational consequences in late-time cosmology, see, e.g., 

[20], mm- 

On the other hand, sources of matter and radiation with typical inhomogeneity 
scales less than rc have somewhat different properties. These are easier to discuss 
for a Schwarzschild source - a spherically-symmetric distribution of matter of the 
mass M and radius ro, such that tm < tq -C Tc (r^ = 2GnM is the Schwarzschild 
radius). For such a source a new scale, that is a combinations of Vc and tm, emerges 

^It is also possible to modify the theory at the linearized level so that the conventional pertur- 
bative expansion is well-behaved [151 [16] , [17] , [18] . 
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(the so-called Vainslitein scal^ [13]: 

r, = {rMrlf''' . (21) 

Below this scale the predictions of the theory are in a good agreement with the GR 
results. Above this scale, however, gravity of a compact object deviates substantially 
from the GR result. Note that is huge for typical astrophysical objects. An 
isolated star of a solar mass would have r^, ~ 100 pc. However, if we draw a sphere 
of a 100 pc radius with the Sun in its center there will be many other starts enclosed 
by that sphere. The matter enclosed by this sphere would have even larger r*. We 
could draw a bigger sphere, but it will enclose more matter which would yield yet 
larger and so on. An isolated object which could potentially be separated from 
a neighboring one by a distance larger than its own is a cluster of galaxies. For 
typical clusters, ~ (/ew Mpc) is just somewhat larger than their size. 

The above arguments suggest that interactions of isolated clusters may be differ- 
ent in the DGP model. On the other hand, at scales beneath a few Mpc or so, there 
will be agreement with the GR results with potentially interesting small deviations. 

For simplicity, we discuss below these issues for a single isolated Schwarzschild 
source. There exist in the literature two different (but both partial) solutions for the 
Schwarzschild problem in the DGP model. The first one is based on approximate 
expansions in the r <C and r ^ regions [21 [131 E] (see also [MJ)- We call this 
set of results the perturbative Schwarzschild (PS) solution. The second one [32] is 
a solution on the brane that interpolates smoothly from r ^ to r ^ Tc ^ r*, 
and is non-analytic in the either parameters used to obtain the PS solution. We call 
this the non-perturbative Schwarzschild (NPS) solution. What is certain, is that 
at observable distances both solutions are in good agreement with the GR results, 
but predict a tiny and potentially measurable deviations from GR [331 [31] (see also 



It is important to understand which of these two solutions, if any, is physically 
viable. Since neither of the two were solved completely in the entire 5D space-time, 
a first step to discriminate between them would be to look closely at the predictions 
that could by tested observationally. This was discussed in detail in Ref. [30]- We 
briefly summarize some of the results. 

Let us start with the Newton potential ^p{r). The result for r ^r^, leads: 



-^^ = -y-^<r^[j) +•••' (22) 

where a = ±v^ and /? = for the perturbative solution (PS), while a = ±0.84 and 
/3 = 3/2 - 2{V3 - 1) ~ 0.04 for the non-perturbative solution (NPS) [32] (two signs 
for a correspond to the tow different branches of solutions, see discussions below). 



similar, but not exactly the same scale was discovered by Vainshtein in massive gravity [12] 
hence the name. 
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The deviation from 4D gravity at r -C gives rise to the additional perihehon 
precession of circular orbits [331 El] (see also [53] for comprehensive studies of these 
and related issues) . In a simplest approximation this effect is quantified by a fraction 
of the deviation of the potential from its Newtonian form 

(23) 

This can be used to evaluate an additional perihelion precession of orbits in the 
Solar system [33l IMj^ . The e ratio is somewhat different for the non-perturbative 
solution (NPS) as compared to the perturbative solution (PS) used in Refs. |33| [3l]. 
This difference has been calculated |36j : 

/ N 0.04 

^ ^ 0.59 ( - ) . (24) 
The perihelion precession per orbit is 

o Q I I / \ 3/2 / N 0.04 

Ay, = 27r + _ _ . (25) 

The second term on the RHS is the Einstein precession, and the last term arises due 
to modification of gravity. For the PS this was first calculated in Refs. [33l [M] ; the 
solution fl25|) is written for the NPS and is somewhat different. 

For the Earth-Moon system r ~ 3.84 x 10^° cm and rf"''*'' ~ 6.59 x 10^^ cm; as 
a result the ratio in fl2^ is approximately 0.48. Therefore, the predictions of the 
non-perturbative solution for the additional perihelion precession of the Moon is a 
factor of two smaller than the predictions of the perturbative solution. The result 
of fl25l) for the additional precession (the last term on the RHS) is =f0.7 x 10~^^ 
(the plus sign for the self-accelerated branch). This is below the current accuracy 
of 2.4 X 10~^^ |37j, but could potentially be probed in the near future [38] . 

A similar calculations can be performed for the anomalous Martian precession 
[51] . For the Sun-Mars system we use rsun-Mars = 2.28 x 10"^^ cm and rf"" = 
4.9 x 10^° cm. The additional precession of the Mars orbit is ~ =f1.3 x 10^^^, which 
should be contrasted with a potential accuracy of the Pathfinder mission ~ 9 x 10^^^. 



4 The quantum gravity scale 

We have discussed in the previous section that the value of the 5D Planck mass M^, 
is restricted by the requirement that rric = 2M^/Mpy be of order Hq, leading to the 
value of the bulk "quantum gravity" scale ~ 100 MeV. 

■^Note that in the leading order of the relativistic expansion the answer is given by the correction 
to the Newtonian potential, while the correction to the rr component of the metric is not important. 
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Here we discuss a slight modification of tlie DGP model which retains the most of 
the important properties of the original theory and yet allows to relax the constrain 
on the bulk quantum gravity scale M^. 

For simplicity we first discuss this for the scaler example. The main idea is that 
the strength of the 5D kinetic term could depend on the y coordinate, so that its 
value is small on the brane but is large off the brane. We could parametrize this as 

dxdyF{x,y){dA(t){x,y))\ (26) 

where the new scalar F (a "dilaton") is assumed to have an x independent profile 
in the y direction such that ^ 0) = mcM|j/4, while 7^ 0) sets the bulk 
quantum gravity scale, which is unrestricted and could be as large as M-p\. 

Let us perform these calculations more carefully. To account for the above prop- 
erties, we introduce an additional term into the action (IT^ which is just a opposite 
sign 5D kinetic term peaked on the brane. To make things tractable, we smear 
the brane, that is, instead of the Dirac function 5{y), we use its regularized version 
5{y) S{y) = 7r-'^e/{y^ + e^), with e ^ 0. The term that we'll be adding to ([12]) 
then reads: 

— j dxdyS{y){dA(P{x,y))\ (27) 

With this term included the variation of the action 6S = with the appropriate 
boundary conditions gives: 

- (Mi, - M^)5{y)dl<P - Mldlct> - dy {{Ml - MH{y))dyct>) = JS{y) . (28) 

Next we take the integral of both sides of this equation w.r.t. y from —e to +e, and 
then turn to the limit 

M ^ 0, e ^ 0, M'^/e ~ M^6{0) = > Ml , (29) 

where we keep M fixed, and its value somewhat lower than M^,. The resulting 
equation reads: 



-dl<t.-''-^^-J^dyC^ = J/Ml,. (30) 



Finally, introducing 



^ 2(M3-M3) 

M|i ' ^^^^ 



where the positive numerical value of nic will be tuned to the Hubble scale today 
rric ^ Hq 10~^^ GeV, we get the desired junction conditions 

-dl<j)-m,dy<j)= J/M^,. (32) 
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Two comments. First, the wrong-sign kinetic term (1271) that is peaked only on the 
brane is dominated by the large positive 4D kinetic term in (IT^ . proportional to 
M|j > M^. Second, the number of adjusted parameters here is the same as in DGP: 
In f|T2l) one should tune the value of M* such that the ratio 2M^/Mpj is of order 
Hq, while in the action (1271) the value of < Mpi can be arbitrary, as long as one 
tunes the value of M so that fl3Tl) is of order Hq. 

In the case of gravity, to which we turn now, similar considerations can applied. 
As before, we smooth out the brane by replacing 6{y) —>■ S{y), and think of the 5D 
EH term to have a profile due to the "dilaton" field 

j d^xdy^5F{x,y)nig5), (33) 

such that the 5D gravitational coupling on the brane is strong, while it becomes 
weak in the bulk. 

The above construction could be parametrized by adding the following boundary 
(worldvolume) term to the DGP action: 

y d'xdy6{y)^,n{g,) = -— j d'x dy5{y)^ N [R + - K^) ,(34) 

where the r.h.s. of (1341) is obtained by using the standard ADM decomposition. The 
total action in the ADM formalism reads: 

Smod = ^ J d^x^,R{g,) + d'xdy^N{R + K^-Kl) 

__y d^xdy5{y)^N{R + K^-K%) , (35) 

where it is implied that the 4D EH term is also smeared over the same scale as the 
5D term0. The equations of motion are straightforward to derive from fl35|) . The 
{/iz/} equation in the bulk, and {/^S} and {55} equations read as follows: 

G'g) = for ||/|>e, (36) 
{Ml - MH{y)) {V^K,, - V,K) = , (37) 
{Ml - M^5{y) ){R-K^ + K^.K^") = . (38) 

As in the scalar case, we will be looking at this theory in the limit (1291) . The above 
equations reduce to ([MH])- 

The Israel junction condition across the brane gets modified because of the new 
term in (l35ll . In the limit (!29|) this condition reads: 

2{M^ — M^) 

- -^^2 (^M- - 9^.vK) = StiGnT^, {x) . (39) 

^For the regularization of 4D and 5D EH terms, see, [39] and [17], respectively. 
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If M = 0, as in ©, we get back the result §^ with = 2M^/M^^. However, M 
does not have to be zero. For an arbitrary value of M^, we tune the value of M so 
that the crossover scale (I3T1) . which appears in ( |39l) . is adjusted to the value of the 
present-day Hubble scale mc ~ -f^o ~ 10~^^ GeV. Hence, fl39|) recovers fH6l) . 

What we have shown is that for M -C Mpi the junction condition is not mod- 
ified as compared to DGP. The 5D gravitational constant in the bulk equations of 
f l55]) would change, though, if we were to consider sources extending into the bulk. 
However, our primary interest is in the sources localized on the brane, and for those, 
the new model (135!) with M <^ Mpi. recovers the results of DGP. 

5 Cosmology 

Let us turn to the cosmological solutions. To this end we consider distributions of 
matter and radiation that are homogeneous at scales ^ Vc- Therefore, the compli- 
cations due to non-linear dynamics outlined in the previous section do not apply to 
these source^. 

The metric is parametrized as follows: 



There are two branches of solutions that are labeled by an integer e = ±1 [T9]: 



P{t,y) = l + e\y\ Q^t,y) = a{t) + e\y\Va^ + k, S(t,y) = l. (41) 



Here we included a nonzero spatial curvature k. With this Ansatz, the Friedmann 
equation on the brane follows from ([3]), and can be expressed in terms of the 4D 
Hubble parameter H = a/a. For a simplest case of a brane without matter/radiation 
on it, and with k = the equation reads [ 19] 



The case with e = +1 admits a dS solution with H = rric. This is called the self- 
accelerated solution, as it gives rise to the accelerated expansion of the universe due 
to modified gravity. We emphasize that the minus sign between the two terms in 
fH2|) is guaranteed by the choice of the positive sign in front of the terms in fHTl) that 
are proportional to If we were to choose the latter signs to be negative, we would 
have obtained the Friedmann equation + mc\H\ = 0, which does not admit the 
dS solution. The latter corresponds to the choice e = — 1, and is referred to as the 
conventional branch. There is no acceleration produced by modified gravity on this 
branch. It has to be obtained, e.g., by introducing a small cosmological constant on 
the brane. An interesting observation concerning this branch is the following [IQ]: 

^Although, those complications will be relevant to perturbations about the background cosmo- 
logical solutions. 



ds^ = -P^{t, y)dt^ + Q'^{t, y)-f,jdx'dx^ + T?{t, y)dy^ . 



(40) 



- em^\H\ = 0. 



(42) 
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the cosmological expansion from a 4D point of view looks as if it had the equation 
of state parameter less than —1. This is because gravity turns into the 5D regime at 
scales H ~ m^, which implies less "deceleration" of matter due to attractive gravity, 
and this appears to a 4D observed as "faster" acceleration |40j . 

It is straightforward to introduce matter/radiation on the brane. Following [19] 
we obtain the Friedmann equation 

which should be amended by the conventional conservation equation for the fluid 
of density p and pressure p: p + 3H{p + p) =0. The latter being a result of the 
matter stress-tensor conservation 'W^T^^ = 0, which can be verified, e.g., by taking 
a covariant derivative of both sides of and using ([7]). 

When H ^ rric, i.e., in the early universe, eq. (H3|) reduces to the conventional 
Friedmann equation. However, at late times, when H ~ rric, the cosmological 
evolution changes. On the conventional branch (e = — 1) it turns into the expansion 
driven by the 5D force law, and on the self-accelerated branch (e = 1) it turns into 
the dS-like expansion. 

How about perturbations on the cosmological solutions? We start with the con- 
ventional branch. Here, the small perturbations about the Minkowski background 
are stable [2j. However, the conventional perturbative expansion breaks down near 
realistic sources, as was discussed in the previous section. As a result, the model 
exhibits the strongly-coupled behavior already at the classical level [13] in the do- 
main where the extrinsic curvature square is > m^. The same applies to the self- 
accelerated branch. Because of this, the question of stability of the self-accelerated 
solution [13] - [16], which happens to be in the strongly coupled regime due to its 
curvature, becomes difficult to address within the perturbative approach [HI [iSj^ . 

Luckily, certain non-linear solutions have been found. This is a case for a 
Schwarzschild-like solution, for which the 4D metric was exactly obtained [32], and 
for the Domain Wall solution for which the full 5D metric was found in Ref. [19]. In 
both cases, the mass (tension) of the solution gets screened by gravitational effects, 
and these sources on the self-accelerated background look as if they had a negative 
net 5D mass (tension)0. This suggests that the self-accelerated background may not 
be problem-free in the full non-linear theory. 

What is a root-cause of this behavior? The self-accelerated solution exists only 
for a certain choice of the sign of the extrinsic curvature, and this choice is such that 

^One could look at this breakdown in terms of the scalar "conformal" mode. In massive gravity 
this mode decouples from the rest of the modes in a certain limit [41]. In the DGP model, however, 
such a decoupling does not take place [H]. 

^In contrast with this, screening of the similar sources on the conventional branch of DGP, 
leaves them with positive 5D mass (tension) [32l |49] . 
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it requires a growing metric in the bulk. For instance, in a simplest spatially-flat 
case, the full 5D metric of the self-accelerated solution takes the form |19j : 

ds^ = (1 + H\y\f {-dt^ + e'^^'dx^} + dy^ , (44) 

where y is the 5th coordinate and H denotes the dS expansion rate of the 4D 
worldvolume (the latter is labeled here by Cartesian coordinates {t,x), and we use 
the 5D coordinate system in which the brane is located at y = 0.). The unusual 
feature of the above metric is that it grows in the bulk, even though the worldvolume 
metric is that of dS space. A linearly growing metric, similar to dS]), would have 
been produced by a negative tension 3-bran^. The growing metric flUl) imprints its 
"negative" effects on the brane worldvolume through the extrinsic curvature, giving 
rise to the solutions mentioned in the previous paragraph. 

Can this problem be cured? It turns out that one can modify the DGP equations 
in such a way that the new system admits a background that is equivalent to the 
self-accelerated solution on the 4D brane, however, differs from it in the bulk. 

The new solution that we will discuss takes the form: 

ds^ = (1 - H\y\f {-dt^ + e^'^'dx^} + dy"" . (45) 

In order to obtain this solution one needs to flip the sign in front of the extrinsic 
curvature term in one of the DGP equations, keeping the rest of the equations 
intact. Below will discuss how such equations can be obtained by modifying the 
DGP action. 

The solution fHS]) is formally identical to that for a 3-brane endowed with a 
positive 4D cosmological constant (brane tension) which is embedded in 5D empty 
space in 5D GR [50j. However, unlike the latter, the worldvolume expansion in the 
present case ( 145|) is due to modified gravity, while the 4D cosmological constant is set 
to zero. This difference is what is responsible for the modified Friedmann equation, 
and distinct cosmological evolution on the self-accelerated background. 

The bulk space in both (jH]) and fH5l) is locally equivalent to 5D Minkowski space. 
In the chosen coordinate system the solution psl) encounters the Rindler horizon at 
\y\ = H~^. However, an analytic continuation beyond this point can be performed 
by employing new coordinates. In that coordinate system the brane (with closed 
spatial sections) can be regarded as a 4D dS bubble that is first contracting and 
then re-expanding in 5D Minkowski space. 

The proposal of Ref. [9] is to add a new term on the brane worldvolume such 
that the sign in front of the second term on the l.h.s. of Q would flip. In other 
words, we introduce a new set of equations in which ([3]) is replaced by: 

G^,y + mc{K^^ - g^uK) + 7^^.^ = 87rG'ArT^^(a;) , (46) 

while all the other equations ([SDjdZD and ([H]) remain intact. Here, in a simplest case 
S^j, = {G^^ - \g^,u{K'^ - Klf^) -2{K^K^a - KK^^)), and the small coefficient 7^1 

^This could be in, e.g., 5D Minkowski or Anti de Sitter (AdS) space-time. 
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is in general nonzero for a regularized brane width. An important property of the 
tensor S^,^ is that it equals to G^^ on the self- accelerated solution given below, while 
it contributes in a non-trivial way to perturbations about it. 

The action functional that gives rise to this new set of equations happens to 
be the one that we have already discussed in the previous section fl35|l . The only 
difference is that in the junction condition fl39l) we need to flip the sign of the 
coefficient in front of the extrinsic curvature terms. This can be achieved if we 
choose M > and set 

_ 2(M^ - M^) 
m. = • (47) 

The metric for the self-accelerated solution of the new system of equations fH6l 
E) reads: 

P{t,y) = l-\y\ Q{t,y) = a{t)-\y\Vd^, S(t,y) = l, (48) 

where we have chosen a negative sign in front of the terms proportional to \y\. Let 
us now see how the change of the positive signs in the metric (HTl) into the negative 
signs in (HHI) changes the value of the extrinsic curvature evaluated at ?/ = 0+. On 
the solution (HSj) . A*"^ = 0, = 1, and K^^ = Hence, at ?/ = 0+ the 

components of the extrinsic curvature tensor evaluated on the solution (1481) equal 
to minus the corresponding components evaluated on fHTl) . Therefore, substitution 
of fHTj) into ([3]) would give the same equation as the substitution of psj) into P6l) . 
The corresponding Friedmann equation on the empty brane, which now follows from 
( H6l) instead of (I3l), reads: 



k /„„ k 



+ ^ = m,\jH^ + ^. (49) 



For = this coincides with ( 1421) and gives the spatially-fiat dS solution with 
H = rric (j45l) . For general k the solutions are: 

ds' = (1 - H\y\f {-dt' + a\t){dx' + Sl{x)dn')] + dy^ , (50) 

where H = rric and k = —1,0,1 corresponds to the open, fiat and closed spatial 
slicing of 4D dS space, for which Sk{x) = sinhx, x, sinx, respectiveljl^. 

The solution fl50l) should satisfy all the bulk equations ([6] -[8]), since in the bulk 
it is locally equivalent to Minkowski space. We checked by direct substitution that 
fl5U]) solves Eqs. (P-[H]) too. The matter /radiation density can also be introduced 
as described above. The Freedman equation coincides with (143!) with e set to 1. 



^There are two other solutions to For fc = one finds the H — flat solution. For fc = — 1 
one finds the Milne solution a{t) = t. 
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It still remains to be shown that the modified model described by (!35|) gets read 
of all the negative mass states that may appear in the self-accelerated solution, this 
question will be discussed in [S3]. I just point out that an additional benefit of 
the new term in fl35l) is that it allows to relax the constraint on the bulk gravity 
scale. The latter can take an arbitrary value below Mpi. This opens a window for a 
possible string theory realization of this model, or its D > 5 counterparts [5], [7] (for 
earlier proposals see 
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